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Foreword 

This report has been prepared to expedite early domestic dissemination of the information generated 
under the contract. The NASA program manager is Dr. M.S. Hirschbein. 

Abstract 

Accomplishments are described for a 3-year program to develop methodology for component-spe- 
cific modeling of aircraft engine hot section components (turbine blades, turbine vanes, and burner 
liners). These accomplishments include: (1) engine thermodynamic and mission models, (2) geome- 
try model generators, (3) remeshing, (4) specialty 3D inelastic structural analysis, (5) computation- 
ally efficient solvers, (6) adaptive solution strategies, (7) engine performance parameters/compo- 
nent response variables decomposition and synthesis, (8) integrated software architecture and 
development, and (9) validation cases for software developed. 
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NOMENCLATURE 


C c = Bound on creep strain gradient 

C_ * Bound on plastic strain gradient 

Maximum allowable sum of 

Lower bound for R^ for possible reswshing 

Maximum allowable upper bound for R^ 

= Bound separating remeshing from re-solving 
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b r = Equivalent heat transfer coefficient for radiation 
to casing Btu 
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K = Metal conductivity, Btu 
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[K] = Stiffness matrix for FEM analysis 
P 3 = Compressor discharge total pressure, psia 


Q/A = Heat flux through material, Btu 

hr-ft z 

{R} = Vector of residual forces in FEM analysis 
R^ = i^ residual force 
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I'HTL) * General Electric proprietary 3D transient beat 
transfer analysis computer program 

TOL ® Tolerance on local integration error 

Tg = Hot side metal temperature, ®F 

* Cold side aietal temperature, °F 

T IIHER * liner temperature , *7 

T3 = Compressor discharge temperature, °F 

T4 s HP turbine rotor inlet temperature, °F 

^41 ~ 2 P turbine inlet gas temperature, # F 

T_ = Temperature at position ij, °F 

T mi1 = Co-lbuator *etal temperature at position 
ij, °F 

t * Material thickness. Ft 
At i = Current time subincrement 
At i+1 = Next time subincrement 
W 41 * Turbine airflow, #/sec 
11X11 = Vector norm 

TC = Cooling effectiveness , dimensionless 

tt » ®» q, r = Temperature dependent material 
creep parameters 

{S} = Vector of displacements is FEU analysis 
(AS e ) #llowable s Maximum allowable stress change 
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(AS ) = Max im u m change in stress occurring in 

e max 

the current time step. 

- Creep strain, m/m 
Ep * Plastic strain, in/in 
E^ - Total strain at point i, in/ in 
E^ = Creep strain at point i, in/ in 
E^ e = Elastic strain at point i, in/ in 
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E^ = Plastic strain at point i, in/ in 

E . 1 = Thermal strain at point i, in/ in 

E* = Second derivative of the inelastic 
strain rate 

(AEg) a j^ow«ble * M a xi mum allowable inelastic 
strain increment. 

(AEg) aax = Maximum inelastic strain increment 
occurring in the current time step 

N c = Nominal cooling effectiveness, dimen- 
sionless 

N cs = Cooling effectiveness at specified span, 
dimensionless 

N = Cooling effectiveness at midspan, dimen- 

COS . * 

sxonless 
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1.0 INTRODUCTION 


Modern jet engine design imposes extremely high loadings and temperatures 
on hot section components. Fuel costs dictate that minimum weight components 
be used wherever possible. In order to satisfy these two criteria, designers 
are turning toward Improved materials and innovative designs. Along with 
these approaches, they also must have more accurate, more economical, and more 
comprehensive analytical methods. 

Numerous analytical methods are available that can, in principle, handle 
any problem that might arise. However, the time and expense required to 
produce acceptable solutions is often excessive. This program addresses this 
problem by setting out a plan to create specialized software packages which 
will provide the necessary answers in an efficient, user-oriented, streamlined 
fashion. Separate component-specific models will be created for burner 
liners, turbine blades, and turbine vanes using fundamental data from many 
technical areas. The methods developed will be simple to execute, but they 
will not be simple in concept. The problem is extremely complex and only by a 
thorough understanding of the details can the important technical approaches 
be extracted. The packaging of these Interdisciplinary approaches into a 
total system must conform to the modular requirements for useful computer 
programs. 

The overall objective of this program was to develop and verify a series 
of interdisciplinary modeling and analysis techniques that have been 
specialized to address three specific hot section components. These 
techniques Incorporate data as well as theoretical methods from many diverse 
areas including cycle and performance analysis, heat transfer analysis, linear 
and nonlinear stress analysis, and mission analysis. Building on the proven 
techniques already available in these fields, the new methods developed 
through this contract have been Integrated to provide an accurate, efficient, 
and unified approach to analyzing combustor burner liners, hollow air-cooled 
turbine blades, and air-cooled turbine vanes. For these components, the 
methods developed predict temperature, deformation, stress, and strain 
histories throughout a complete flight mission. 
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This program, to a great extent, drew on prior experience. This base of 
experience was invaluable for understanding the highly complex intersections 
among all the different technical disciplines as well as for estimating the 
importance of different engine parameters. In particular, there are four 
specific areas in which experience was especially beneficial. 

First, with the recent Increases in fuel costs, greater emphasis has been 
placed on more accurate solutions for stresses and strains in order to 
understand and improve the durability and life of hot section components; 
conventional linear elastic analyses are no longer sufficient; instead, they 
now provide the boundary values for more refined creep and plasticity 
calculations. These nonlinear analyses are now performed routinely as part of 
the design process at General Electric. This extensive experience with these 
plasticity and creep methods contributed directly to developing component 
specific models. 

Second, advances in 3-D modeling capability are being achieved by the 
concepts developed under the NASA-supported ESMOSS program. ESMOSS concepts 
provided the basis to develop an efficient modeling system for geometric and 
discretized models of engine components. 

Third, the NASA-funded Burner Liner Thermal /Structural Load Modeling 
Program contributed strong support to this program. The specific area 
addressed the transfer of data from a 3-D heat transfer analysis model to a 
3-D stress analysis model and provided the background and framework for the 
data interpolation required for all thermomechanical models in this contract. 

Fourth, over the past 10 years General Electric has developed Internally a 
family of computer programs: LASTS, OPSEV, and HOTSAM. These programs all 

have the common thread of using selected points from cycle data, heat transfer 
and stress analyses, and a decomposition/synthesis approach to produce 
accurate values of temperature, stress, and strain throughout a mission. 

These programs were totally consistent with the overall objectives of this 
program and represented a proven technology base upon which the component 
specific models were developed. Significant advances were made in the 
inclusion of nonlinear effects and the introduction of improved modeling and 
data transfer techniques. 
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The program was organized into nine tasks which were logically separated 
into two broadly parallel activities (Figure 1). On the right of Figure 1 is 
the Component Specific Thermomechanical Load Mission Modeling path. Along 
this path a Decomposition/Synthesis approach was taken. In broad terms, 
methods were developed to generate approximate numerical models for the engine 
cycle and the aerodynamic and heat transfer analyses needed to provide the 
Input conditions for hot parts stress and life analysis. 

The left path, Component Specific Structural Modeling, provided the tools 
to develop and analyze finite element nonlinear stress analysis models of 
combustor liners, turbine blades, and vanes. These two paths are shown in 
more detail in Figures 2 and 3. 

Software Development, Task IV, consisted of planning and writing the 
computer programs for both paths, with the necessary interconnections, using a 
structured, top down approach. 

In the Thermomechanical Load Mission Modeling portion of the program 
( * ask III), a Thermodynamic Engine Model was developed which generates the 
engine Internal flow variables for any point on the operating mission. 

Figure 2 diagrams the tasks. The method for doing this is described below. 

Task V was developing techniques to decompose flight missions into 
characteristic mission segments. In Task VII a Thermomechanical Mission Model 
was developed. This uses the flow variables from the Thermodynamic Model to 
determine metal temperature and pressure distributions for a representative 
combustor liner and turbine blade and vane. 

Individual tasks for the Structural Modeling activity are shown in 
Figure 3. The requirements of Software Design, Task II, were factored into 
Task VI, the evaluation of the structural analysis methods which were selected 
for evaluation in Task I. Task VIII provides the capability for structurally 
modeling current state-of-the-art combustor liners and hollow turbine blades 
and vans, given the defining dimensional parameters. These parameters were 
chosen to facilitate parametric studies. 
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The component specific models were developed In two steps. In the first 

a geometric model is defined, in the aDDliratinn +h r 
M . in ine a PP Mcation of the Component Specific 

HI s I"V7?7 " ata arC the " trans,errei t0 «• Thermomechanical Load 

temperatures ° Pr ° V,d ® ^ 9e0metry f0r <letern, ' n1 "9 component pressures and 
temperatures. Thus, a data transfer link was developed to do this In Task IV 

r :::;”:::’ 0 ;:::- r capab,my f ° r *-^,1 ’ 

discretized, finite element model was a part of Task VIII. At this point 

;; t ,,nk “* *" « -eeded to transfer the complnt 

temperatures and pressures from the Thermomechanical Load Model to the finite 
e ement model, Interpolating the data as needed to define nodal temperatures 
nd pressures. This was completed in Task IV. 
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Figure 3. Component Specific Structural Modeling. 
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2.0 TECHNICAL PROGRESS 


2.1 TASK I - LITERATURE SURVEY 

The first task of this program was to perform a literature survey of 
available methods, techniques, and solution strategies that could be used to 
geometrically model, display, and structurally analyze burner liner, turbine 
blades and vanes. NTIS, NASA, OTIC, and internal General Electric Company 
documents were searched. As a result of this survey, 85 papers and 8 books 
and procedures were discovered with pertinent information. As a result of 
evaluating this information, recommendations were made on the technology to be 
incorporated into the program and approved by the NASA Program Manager. 

2-2 TASK II - DESIGN OF STRUCTURAL ANALYSIS SOFTWARE ARCHITECTURE 


The software architecture was designed using the methodology developed 
for the ESMOSS program. This development was carried out by a team whose 
members provided expertise in all of the pertinent areas. The architecture 
approved by the NASA Program Manager is contained in Appendix A. In addition 
to the program architecture, the preprocessor and postprocessor attributes are 
defined. 

2.3 TASK III - THERMODYNAMIC ENGINE MODEL 

The Thermodynamic Engine Model has been developed as a simple calculation 
tool which will take as inputs the three variables altitude (h), Mach number 
(M), and power level (PL) over the allowed flight map of an engine, as shown 
in Figure 4. In addition, ambient temperature deviations from the standard 
atmosphere, airframe bleed air requirements, and engine deterioration can also 
be included as part of the input to the thermodynamic model. For each input 
condition specified by h, M, and PL, the thermodynamic model will calculate 
gas weight flow (W), temperature (t), and pressure (p) at selected aerodynamic 
engine stations as needed to determine component thermal loadings. These 
stations are shown in Figure 5. 
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The technique for developing a thermodynamic engine model is shown in 
Figures 6 and 7. The engine to be analyzed must be defined thermodynamically 
by an engine cycle deck (computer program) that can be run to generate the 
internal flow variables at the chosen aerodynamic stations (Figure 6). To 
encompass the complete engine operating map (Figure 4), 148 operating points 
are chosen and w, t, and p are calculated using the cycle deck for the 
selected stations as well as and N 2> the fan and core speeds. From 
this station data, an engine performance cycle map is constructed. This is 
essentially a set of three-dimensional data arrays that map the station data 
(w, t, p, N 1 , and N 2 ) onto the engine operating map (Figure 4). Given an 
arbitrary operating point defined by h, M, and PL, it is then possible in 
principle to interpolate on the engine performance cycle map to determine 
station data. In practice the station parameters are nonlinear functions of 
the input parameters, and considerable effort is needed to develop these 
multidimensional interpolations. The computer programs used to generate the 
engine performance cycle map from the engine cycle deck output has been 
developed as part of Task III. The functioning of the thermodynamic engine 
model is shown in Figure 7. Given an engine mission, as shown schematically 
in Figure 8, it can be defined by values of the input variables h, M, and PL 
at selected times through the mission. Using these input variables and 
the engine performance cycle map, an interpolation program developed in 
this effort calculates engine station parameters throughout the mission 
(Figure 7). These are then used to define station mission profiles of w, t, 
p, N 1 , and N 2 as functions of time at each aerodynamic station. These 
station mission profiles then become the input to the thermomechanical engine 
model . 


2.3.1 Detailed Specification and Requirements 

The first step in this task was to develop detailed specifications and 
requirements for the thermodynamic engine model software. This specification, 
which was approved by the NASA Program Manager, is presented in Appendix C. 


Mach No 
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Power Level 

(GI * Ground Idle) 


Figure 4. Engine Operating Map. 
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2.3.2 Model Design and Development 


Based on the detailed specifications, the thermodynamic engine model 
software was generated and a set of 148 performance cases was obtained to load 
this model. The next task was to establish interparameter interpolation 
functions. To assist in this effort, 103 special cases that would maximize 
Interpolation errors were chosen from the cycle deck. 

A master Interparameter linearity study was executed to evaluate 
interpolation functions. A computer program (STATPAC) was available that 
could take 30 input performance parameters, perform transformations on the 
data, and generate crossplots of the transformed function. The linearity of 
these crossplots was the criterion of excellence in the selection of 
interpolation functions. One hundred validation cases were run with 30 
parameters each, giving 3000 individual comparisons. Sixty-three additional 
performance cases were used to perform the interparameter linearity study for 
the Mach number and altitude control variables. 

Based on the above program, a set of Interparameter interpolation 
functions and transformation functions were defined and encoded in the TOE 
model software. The accuracy that can be achieved with these is excellent. 

As a final "trial run," this model was tested against the CFM56 engine flight 
conditions. A TOE User's Manual has been written. 

2.4 TASK IV - SOFTWARE DEVELOPMENT 

2.4.1 Software System Overview 

The Component Specific Modeling Program is composed of several software 
packages. These primary constituent software packages are the Thermodynamic 
Engine Model, the Thermomechanical Loads Model, the Component Specific 
Geometric Models, and the FEM Structural Analysis Code. Each of these 
software packages exist and function as separate Independent codes. It has 
been the primary objective of the COSMOS program to concentrate the effort of 
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the diverse technologies Into a single software system. Reference to the 
system design of Figure 9 provides the reader a good overview as to the 
structure of the COSMOS software system. 

In brief, a main executive controls overall program function and data 
flow, while separate functional modules perform defined tasks. 

2.4.2 Thermodynamic Engine Model 

The Thermodynamic Engine Model has been completed and installed at NASA 
Lewis. The model has been developed as a simple calculations tool which will 
take as inputs the three variables: altitude (h), Mach number (M), and power 

level (PL) on the allowed flight map of an engine, as shown in Figure 5. In 
addition, ambient temperature deviations from the standard atmosphere, 
airframe bleed air requirements, and engine deterioration can also be included 
as part of the input to the thermodynamic model. For each input condition 
specified by h, M, and PL, the thermodynamic model will calculate gas weight 
flow (w) , temperature (t), and pressure (p) at selected aerodynamic engine 
stations as needed to determine component thermal loadings. These stations 
are shown in Figure 5. 

The technique for developing a thermodynamic engine model is shown in 
Figures 6 and 7. The engine to be analyzed must be defined thermodynamically 
by an engine cycle deck (computer program) that can be run to generate the 
internal flow variables at the chosen aerodynamic stations (Figure 6). To 
encompass the complete engine operating map (Figure 4, 148 operating points 
are chosen and w, t, and p are calculated using the cycle deck for the 
selected stations as well as and N 2 , the fan and core speeds. From 
this station data, an engine performance cycle map is constructed. This is 
essentially a set of three-dimensional data arrays that map the station data 
(w, t, p, N.j , and N 2 ) onto the engine operating map (Figure 4). Given an 
arbitrary operating point defined by h, M, and PL, It Is then possible In 
principle to Interpolate on the engine performance cycle map to determine 
station data. In practice the station parameters are nonlinear functions of 
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the input parameters, and considerable effort is needed to develop these 
multidimensional interpolations. The computer programs used to generate the 
engine performance cycle map from the engine cycle desk output has been 
developed as part of Task III. The functioning of the thermodynamic engine 
model is shown in Figure 7. Given an engine mission, as shown schematically 
in Figure 8, it can be defined by values of the input variables h, M, and PL 
at selected times through the mission. Using these input variables and the 
engine performance cycle map, an Interpolation program developed in this 
effort will calculate engine station parameters throughout the mission 
(Figure 7). These are then used to define station mission profiles of w, t, 
p, N.| , and as functions of time at each aerodynamic station. These 
station mission profiles then become the input to the thermomechanical engine 
model . 

2.4.3 Thermomechanical Loads Model 

The Thermomechanical Loads Model is a computer model based on various 
types of correlations previously developed within GE. For our base program, 
work has been completed on a generalized procedure that has been established 
to predict surface temperatures and pressure loads for a rolled ring 
combustor. This procedure has been developed using all available data arising 
from both experimental test and analytical methods. 

Input for such a model is relatively simple and limited in scope. 
Geometric information, coupled with heat transfer coefficients and pressure 
scaling constants, is all combined to define the resulting metal temperatures 
and pressure loads. A cross-mesh matching algorithm is then employed to map 
the resulting temperatures and loads to the FEM model generated by the 
component specific geometric software. 
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2.4.4 Structural Analysis Code 


2. 4. 4.1 Overview of Attributes & Capability 

The structural analysis code developed under this program has been 
designed to solve complex high temperature thermomechanical problems involving 
cyclic creep and plasticity. The development of this program has been 
directed toward its application to the complex geometries, loads, and high 
temperatures associated with the components within the hot section of gas 
turbine engines. The finite element code was designed and developed for this 
contract and is somewhat unique in both its structure and content 
(Figure 10) . 

The analysis code is a combination of what we have deemed the most 
desirable features of many analysis codes into a single entity. Some of these 
features are presented here. 

2. 4. 4. 2 Preprocessing 

Preprocessing has been made an integral part of the analysis code. With 
such capability, the user can instruct the program to read all input and check 
for errors without having to execute a separate preprocessing type program. 

The preprocessing function essentially follows the data flow through the 
analysis code and provides the user a data summary and, in turn, alerts the 
user to any detected analysis difficulties. 

2. 4. 4. 3 Bandwidth Reduction/Optimization 

An automated bandwidth reduction algorithm (Gibbs/Poole/Stockmeyer) has 
been incorporated into the analysis code. This feature may be used as either 
a selected or a default option. Since most FEM solution algorithms are keyed 
to bandwidth, a reduction thereon provides the user a more efficient and cost 
effective solution. The bandwidth reduction algorithm plays an important role 
in the ability to efficiently and effectively Implement remeshing and mesh 
refinement techniques/algorithms. 
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2. 4. 4. 4 Internal Numbering System 


An internal node and element numbering system is contained within the 
analysis code. Such a system is extremely useful and affords both the user 
and the programmer great utility. The user is now free to use numbering 
systems of convenience and allows for the removal or editing of nodes and 
elements in the model without affecting the code's overall operation or 
efficiency. 


2.4.4. 5 Internal Data Transfer 

Much effort was directed towards the efficient storage and retrieval of 
data within the analysis code. Working information has been grouped into data 
packets to minimize file manipulation and processing time. 

2. 4. 4. 6 Solution Techniques 

For reasons of efficiency and problem application dependency, provisions 
for alternate solution techniques have been included within the analysis 
code. Currently contained in the analysis code is a COLSOL algorithm with 
provisions therein for SESOL and FRONTAL techniques. 

2. 4. 4. 7 Constitutive Models 

The simplified constitutive model and Haisler & Allen's are used to 
represent inelastic material behavior. Provisions are included in the code to 
incorporate other constitutive models as they are developed or as they need to 
be included. 
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2. 4. 4. 8 Element l_ihrary 


The twenty noded isoparametric is the element of choice for the 
demonstration code. The performance of the element for elastic analysis as 
well as inelastic material applications is well documented. The element 
stiffness matrix can be computed using either second or third order Gaussian 
quadrature, but stress recovery is available only at the second order 
integration points. The eight and sixteen noded isoparametric elements will 
be incorporated into the code during the follow-on program. Their inclusion 

will be direct and straightforward owing to the modularity of the overall 
code. 


2-4. 4. 9 Postprocessing 

All postprocessing type functions excluding remeshing and mesh refinement 
og c are provided via the commercially available, general purpose PATRAN 

program. All data files required by the PATRAN program are created within the- 
analysis code and are ready for input. 
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^•5 TASK V - MISSION MODEL DEVELOPMENT 


2 * 5,1 -- P ° nent Temperature an d Pres sure Pecomnn<iti on an ri Synthesis 

Based on the developments described below, a Component Temperature and 

Manaaer 6 ^° mP ° Siti ° n 3nd S y" thesi s Plan was approved by the NASA Program 
Manager. This plan is outlined in Appendix B. 

2 ‘ 5 * 1 ' 1 ~ Pustor Liner Temp eratu r e and Press.i™ Decomnns-m™ 
and Synthesis 

An expression for the temperature gradient through the material 
thickness can he derived from the cooling effectiveness, the compressor 

rir,:-?"- ^ the «« temperature. The temperature 

gradient through the material can be calculated from 


O) 


t h - T r = (mi 

H C K 


where 


^ = Hot side metal temperature, °F 

Tg * Cold side metal temperature, °F 

0/A « Heat flux through material — * u 

hr- ft 2 

t = Material thickness, ft. 


Btu 

ft 2 -hr-° F/ft 


K = Metal conductivity 



The heat flux can be calculated from 

Q/A = ( h c + h r>( T C - t 3 ) 

or It is proportional to (T, . .t i 

Liner V 

0/A = (h + h )(T, , - t ) 

c r Liner 3 ' 


where 


h c = Convection heat transfer coefficient Btu/hr-ft 2 -°F 
h r = Equivalent heat transfer coefficient for radiation 
Btu/hr-ft 2 -°F 

^Liner = liner temperature, °F 


gives 


Substituting the heat flux expression into the 


gradient 


T H ~ T c 


(T Liner “ V 


0» c + hr)t 


K 


using the equation for cooling effectiveness, 

T/. ’ T t ,. 


V = —A Line r 

‘ T *- T 3 


Sn equation for (T Liner - T } ) can be written as follows: 
(T L1ner “ T 3 ) = (1 ' V c ) T 4 + <y c - D T 


( 2 ) 


( 3 ) 


to casing 


equation (1) 


( 4 ) 


( 5 ) 


( 6 ) 
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Substituting Equation 6 into the expression (4) gives 

t H " T 3 0 * c + br)t 

a - y c ) t 4 + (y c - 1) t 3 a i (7) 

The convection term, h c , varies with pressure and thus the gradient through 
the material thickness should be correlated with pressure. 

A THTD analysis was done at several pressure conditions and the 
calculated temperature gradients were plotted vs P3 for several axial 
locations and the results are shown in Figure 11. The locations and 
coordinates are shown in Figure 12. As shown in the figure, the gradient data 
are correlated with pressure. The constants m and b in the equation 


T - 
n 

^Liner 




a - y c > T y * ch c - i) t 3 


= mP 3 + b 


( 8 ) 


are tabulated in Table I. 

Table I. Linear Fit Constants for Equation (8). 
| Location X, inches M b 


1 

0.094 

12.3 x lO- 5 

0.100 

2 

0.438 

14.1 x 10-5 

0.061 

3 

0.654 

9.0 x 10-5 

0.061 

4 

0.854 

10.7 x 10-5 

0.092 

5 

1.114 

28.1 x 10-5 

0.168 
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The spanwise distribution of overall local cooling effectiveness has been 
completed for two different Stage 1 HP turbine blades. The results are 
compared In the attached Figure 13 which Is based on Transient Heat Transfer 
Analysis runs for 15, 50, and 77% of blade span. The ratio of 

n c local /n c 50% unity by d ®fi n ttion at the 50% span location. At 
the other two spans, the h ratio is identified for each of the sixteen 
points around the airfoil. The curves have been terminated at the locus or 
the average h c ratio for each span. 

For the first blade, this procedure defined a single curve for the 
pressure and suction surfaces. However, the other blade is better represented 
by a two-branch curve at the 77% span (Figure 14). We have reviewed the two 
blade designs for possible explanations of this characteristic. There is no 
obvious single cause. It is undoubtedly the combined result of configuration, 
coolant circuitry, the application of film cooling and variations In gas-side 
heat-transfer coefficients. 

It appeared best to allow for incorporating separate curves for the 
pressure and suction surfaces, with freedom to input these curves for 
different blade designs. This is probably the thing to do for the second 
blade defining a separate curve for the pressure surface between Points 2 
and 6 (Fig. 14). Points 7 and 8 appear to be represented quite well by the 
curve for the suction surface. Using the suction surface curve for Points 2 
through 6 could overpredict the temperatures by about 135°F at the 77% 
span. 

2.6 TASK VI - STRUCTURAL ANALYSIS METHODS EVALUATION 

This task was completed during 1983. The selections approved in 2.1 are 
the result of this effort. 
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Figure 14. HPT Blade “B” Cooling Effectiveness. 
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2.7 TASK VII - THERMODYNAMIC LOADS HOPE! 


Thermodynamic Loads Features for TOE Model: 
I. BLADES AND VANES 
A. Input is : 


1. A series of midspan station temperatures at a specified SS 
reference case condition, and for each station a percent cord 
envelope dimension 


2. A set of output percent cord envelope dimensions 


3. A table of percent radial span versus cooling effectiveness 
factors. RF, where: 


RF = 



n 

'ass 


n 

] 

n 

i 


n = cooling effectiveness at specified span dim. 

'cs 

^csm = cooling effectiveness at midspan 

The point density in this table will be such that spanwise linear 
interpolation will suffice. 

B. Calculations 


! 


n 


1. For each temperature a nominal cooling effectiveness is 
calculated from: 


n 


c 


T41-Tm 


T41-T3 
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2. At each new condition, all station n c values are modified by 
a factor, CF, as follows: 

• n c ) 

CF = (1 - n c )REF = ^ ^ T3*T4lj^j.) / (T3^£j.*T4l) )a, 

where a is input or default value. 

3. Output n c values are linearly interpolated based on percent 
cord envelope at adjacent input stations. 

4. For each output radial distance (specified in the input) all 
station temperatures are calculated as follows: 

Tlj = T3 + (1 - n )*CF(T41 -T3)*RF 
C R 


Where i = station index 

j = radial station index 


II. COMBUSTOR 


A. Inputs are: 

1. Metal bulk node temperatures at axial stations on inner and outer 
liners, and for each station, x and y dimensions of node centers 
at a specified reference case condition. Both hot streak and 
average metal temperature values are input. 

2. It is assumed that output node dimensions will match the input 
nodes. 
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3. A set of linear equation constants are input for a sparse set of 
axial locations, identified by "AT station" numbers. Output 

AT locations are same as input locations. 

4. A set of AP scaling constants are input at a sparse set of 
locations, identified by AP station numbers. Output AP 
stations are same as input. 


B. Calculations 

1. A nominal cooling effectiveness is calculated for each input 
combustor average and hot streak temperature. 

2. At each flight condition, each metal temperature is recalculated 
from: 

Tm ij = T3 j +0 " H c )(T41 j-T3j) 

where 1 = station Index and j = flight phase index 

3. At each flight condition, each AT Is recalculated from: 

AT 1j = [T3 j + 0 " n c _)(T41-T3)][b i +m i *P3 ;j ] 

4. At each flight condition, each AP is recalculated from: 

"y * p V K i* (b 1 )’ * T3 j 
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The accuracy of the thermodynamic engine model has been evaluated, 
relative to the engine steady-state performance computer model (cycle deck), 
in a model based on interpolation methods, the maximum error must occur in 
regions where the "distance" from known data points is greatest. Figure 15 
shows the altitude versus Mach number map of stored data points used in the 
thermodynamic engine model (TEM). The X symbols on this map Indicate the 
worst-case points selected for the error evaluation. 

The interpolation logic varies from quadrant to quadrant within a 
rectangle bounded by stored data points. The nearest point is always used as 
the base from which the interpolation process is started, for example. 

Figure 16 shows a typical set of four stored data points. The shaded area is 
a quadrant. The evaluation test point is at the center. The hypothetical 

error magnitude curves drawn from two corners to the center illustrate two 
facts : 

• In any quadrant, the error surface is approximately parabolic in 
shape and maximum at the center of the stored data. 

# The four quadrants have different error surfaces, and a 
discontinuity occurs where they meet. 

Since there can be four different maximum error values at each test 
point, the error analysis was performed four times at each point. The results 
were summarized by a computer program. The right-hand four columns of 
Table II show the error values that exceed the target value which Is listed in 
the center of the heading. The left-hand five columns identify the test 
points. The middle four columns show the accuracy level available before the 
improvements developed in this program were incorporated. The average of the 
absolute errors, the max error, and the number of "Exceedances" are given at 
the bottom of each column. Note that speed, pressure, and horsepower errors 
are expressed as percent of the rated standard day, sea level value. It 
seemed more meaningful to express temperature errors in degrees. 
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Table II. Error Summary. 
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Table III shows a brief summary of the accuracy level achieved. Column 2 
shows the average of all test point errors. Columns 3, 4, and 5 show the 
value that is exceeded 2, 4, and 11 times (IX, 2X, and 5% of the 220 error 
values). Note that all data in this figure refers to the worst-case test 
points. Since the error surfaces are approximately parabolic in shape, the 
average error in each quadrant is approximately half of the maximum error, and 
the overall error is approximately half of the average error listed. 



Table III. 

Validation Case 
Error Analysis. 




Error 

Exceeded 

N Times 


Average 

N = 2 

N = 4 

N = 11 

P 2 

0.03X 

<1X 

<1X 

<1% 

P 3 

0.23X 

1.3X 

1.2% 

1.0% 

FNIN1 

0.49X 

1 .8% 

1.8% 

1.8% 

XN25 

0.17X 

0.7X 

0.6% 

0.5% 

T 2 

0.08% 

A 

tn 

o 

<5° 

<5° 

T 3 

2.1° 

<10° 

<10° 

A 

O 

O 

T 41 

11° 

47° 

47° 

35° 


The maximum error of 47°F listed for T41 may seem large. This was one of 
the most difficult parameters to fit. However, note that 47° is only 2X of 
the rated T41 value (expressed in °F) and the true average error is only 
approximately 0.23X. A User's Manual has been written. 
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2.8 TASK VIII - COMPONENT SPECIFIC MODEL DEVELOPMENT 


2.8.1 Geometric Modeling 

2. 8. 1.1 Combustor Liner Nugget Model 

The objective of this subtask was to develop and Implement an integrated 
procedure that would allow the user to create a 3D finite element model of a 
sector of a combustor liner. The software developed allows the user to 
interactively create the nugget geometry, maps the thermodynamic engine loads 
(temperature and pressure) into the geometry, generates both 2D and 3D finite 
element models and interprets the engine temperatures into the 3D model using 
the average thermal values and the hot streak values. It also automatically 
applies the appropriate combustor pressure to the 3D model and generates the 
symmetric boundary constraints. Finally, it writes this Information in a form 
that is acceptable to the 3D analysis code being used. 

The basic nugget geometry is generated using a "geometry recipe" 
approach. This technique allows the user to generate and control the geometry 
through a small set of geometrical parameters. These parameters include data 
which describe the thickness of the liner, the radii of the nugget and the 
slopes and lengths of key segments in the nugget. Figure 17 shows these 
parameters and their default values for the liner. Once these parameters are 
set by the user, the geometry module generates a 2D finite element model 
representing the liner and the nugget geometry specified. Figure 18 
illustrates a typical 2D model with element numbers. The user is then 
prompted for the thermal/mechanical load input file. This file contains the 
temperature and pressure data produced by the thermal /mechanical load module 
for both the average quantities and the hot streak values. Figure 19 shows an 
example listing of this file. The user is also prompted for the number of 
fuel nozzles and the number of circumferential elements desired. The module 
automatically generates a 3D sector model representing an Included angle equal 
to 180° divided by the number of nozzles. The circumferential elemental 
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Figure 17. Combustor Liner Parameters. 
















X 

_Y_ 

0.034 

. 030 

0.094 

0.030 

0.159 

0.030 

0.208 

0.030 

0.254 

0.030 

0.298 

0.030 

0.340 

0.030 

0.404 

0.030 

0.438 

0.030 

0.584 

0.040 

0.654 

0.040 

0.764 

0.040 

0.854 

0.040 

0.934 

0.044 

1.040 

0.044 

1.074 

0.040 

1.114 

0.030 

1.150 

0.030 


T„ 

T 

H 

A 

520.182 

520.182 

521.284 

. 521.284 

522.834 

522.834 

523.901 

523.901 

524.487 

524.487 

524.935 

524.935 

525.658 

525.520 

526.519 

526.002 

528.827 

526.967 

520.446. 

527.690 

531.376 

528.276 

531.686 

528.758 

530.928 

527.690 

529.826 

527.311 

528.586 

526.967 

527.828.. 

526.726 

527.656 

526.485 

527.518 

526.381 



_^A 

AP 

0.161 

0.161 

0. 

0.263 

0.263 

0. 

0.391 

0.391 

0. 

0.463 

0.463 

0. 

0.485 

0.485 

0. 

0.499 

0.499 

0. 

0.516 

0.505 

0. 

0.523 

0.489 

0. 

0.639 

0.521 

0. 

0.735 . 

0.562 

0. 

0.790 

0.597 

0 . 

1.032 

0.800 

0 . 

1.144 

0.841 

0 . 

1.310 

1.014 

0 . . 

1.481 

1.238 

0 . 

1.461 

1.284 

0 . 

1.541 

1.340 

0 . 

1.613 

1.405 

0 . 


Figure 19. Outputs from Combustor Thermodynamic Loads Model. 
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Figure 22. Typical Program Run. 


2. 8. 1.2 Geometrical Modeling of Turbine Blade and Vanes 


The recipe approach was used very successfully for the combustor liner; 
however, It cannot be used for the blade and vanes due to the significant 
complexity of the airfoil geometry. A typical multiple cavity turbine blade 
is shown in Figure 23. There is no reasonable number of physical parameters 
which could be used to easily describe this geometry. The outside contour of 
the airfoil is a complex curve defined by aerodynamic requirements. In many 
cases, the aero analysis is performed prior to or in parallel with the 
mechanical design, and one of the outputs of the aero analysis is a file of 
points defining the outer contour of the blade at a number of span locations. 

The models developed in this subtask use as input a coordinate file which 
defines the outer airfoil contour for a maximum of 30 spanwise sections. The 
number of points per section can vary with a maximum of 60. These limits are 
parameterized in the code and can be easily adjusted. Figure 24 illustrates a 
series of points for two spanwise sections. 

The code allows the user to specify the number of elements in the spanwise and 
chordwise directions and also allows for multiple elements through the 
thickness of the airfoil. The program will automatically evenly space the 
elements in all three directions or the user can control the biasing of the 
elements. The basic section input can be modified using built in scaling and 
offset logic, and the leading and trailing edge thicknesses can be altered 
through user prompts. The number, size and location of the internal cavities 
can also be controlled by the user. Figure 25 shows an example of a simple 
execution of the code. 

Figure 26 shows a 30 model comprised of 20-noded isoparametric elements of a 
solid airfoil created by this module. Figures 27 and 28 represent the same 
airfoil with two cavities near the leading edge. The plot in Figure 28 is a 
"free edge" plot that clearly shows the cavity as it passes through the 
airfoil. Figures 29 and 30 show a model created with three large cavities and 
three elements across each of the spans. Note that wall thickness of the 
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airfoil was varied along the span for this example. The final example shown 
In Figures 31 and 33 Illustrates how the spanwlse and chordwlse element 
biasing can be used In conjunction with a variable layer thickness. Once 
again, the output of this module Is a file containing 20-noded solid elements 
compatible with the 30 nonlinear code developed as a part of this contract. 
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Figure 24. Points Describing Outer Airfoil Contour. 
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TITLE- 


DATE - 10/2-1 /B6 


TIME » 12.02 


PA3E* 1 


COMPONENT SPECIFIC AIRCIL GENERATOR 

*#*» ***** 0 0000000000 *, 000000 0t.0 00***0*0*0**0000000 

INPUT AIRFOIL FILENAME 
-AER02 

HOW MANY PAIRS OF POINTS IN EACH CROSS SECTION? 

-17 

HOW MANY DIGITIZED SECTIONS 
*2 

TYPE 1 TO MODIFY THE INPUT DATA 

*0 

HOW MANY ELEMENTS ALONG CHORD ? 

» 10 

TYPE 1 FOR USER DEFINED SPACING 

-0 

SECTION 1 UNWRAPPED CHORD 1.8337 
SECTION 2 UNWRAPPED CHORD 1.8316 
HOW MANY ELEMENTS ALONG SPAN ? 

-5 

HOW MANY ELEMENTS ACROSS THE THICKNESS 

•3 

TYPE O FOR THE EQUALLY SPACED METHOD 

TYPE 1 FOR THE MULITPLE COATING (SKIN) METHOD 

•0 

ENTER THE NUMBER OF CAVITY LAYERS 

« 1 

ENTER THE NUMBER OF CAVITY ELEMENT ROWS 
AND THE ELEMENT CAVITY ROW NUMBERS 
•2 3 5 

YOUR UIF IS ON FILE 31 

NO. OF NODES , LAST NODE • 922 922 

NO. OF ELEMENTS - 140 


Figure 25. Simple Example of Airfoil Generator Execution. 
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2-8.2 Remeshing and Hesh Refinpmont 


This area and the next. Self-Adaptive Solution strategies, touch on each 

°! S ? a " y ' S °“ 9ht th,S pr0 9 ram 11 the hes, combination 

both. This involves two major areas of Investigation: the method to be 

use o refine, upgrade, and rearrange the mesh, and the criteria to be used 
to activate this process. 


8 n “ mber ° f WSyS t0 r6f,ne a “»* a better answer: (1) 

one way Is to progressively subdivide a coarse mesh, always retaining all 

prev ous meshes within the finer mesh; (2) , second family of techniques 

tota ly realigns the mesh based on some criteria such as strain energy 

density; (3) a third method Is to leave the mesh unchanged but upgrade the 
order of the elements. 


The first method, progressive subdivision, has certain theoretical and 
computational advantages. If the finite element Interpolating functions used 
meet the requirements for completeness and continuity, convergence Is 
mathematically guaranteed when we refine the mesh by progressive subdivision 
he computational process of remeshing by progressive subdivision Is 
straightforward; however. It guarantees a larger problem to solve. 


For a solution of the finite element system of equations: 
[K] {«} = {F} 


suppose there Is a numerical solution for the displacement, { 6 *}. Then 
the equilibrium or residual force vector is generated: 


{R> = {F} - [K] {w*} 
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A perfect solution would result in this vector containing all zeros. 
Given the finite numerical accuracy of the computer, this is impossible. 
Therefore, a measure of the numerical “goodness" of the solution is to be 
found in how much this vector deviates from zero. Decisions on whether to 
re-solve or redefine the problem can be based on the total and local 
deviations from zero. If a few local degrees of freedom are out of 
equilibrium, this might suggest a local remeshing. If the total equilibrium 
is deficient, this will require remeshing and/or re-solving with greater 
numerical accuracy. 


The decision tree for this is as follows 


1. i If ZR. < C D 

l K 


and all R. < C 


1 


'li- 


the solution is good 


2 . 


If !R i > C R 

and [Number of nodes with C 
then re-solve 


<1L 


< Rj 


< C,. ) > c. 


'1U 


3. If ZR i > C R 

and [Number of nodes with R^ > C R ] < C 
then remesh and re-solve lu 


4. If R i < C R 

but some R. > C_ 

* K iu 

then remesh and re-solve 


where: 

R., = 1 th residual-free vector 

C D - Maximum allowable sum of R. 

R iL 1 

C D = Lower bound for R1 for possible remeshing 
K iu 

C R = Maximum allowable upper bound for R^ 
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°r: an acceptable ‘W'—t solution has been reached. proceed to the 
rir; ^ e ' iStiC ,eVe ’' StraSSeS a " d “"ins are linearly 

heck d Th , y 7 theS<! tW ° " ae<IS l ° b6 evaU,atei1 - * 

of an ei.stT * Stra,n “ MCh Ca,cu,4t, »" <»'"» in an element Is made up 
Of an elastic strain and a thermal strain: 

-e? ♦ «! 

Ill 


gradients for diff ! T"™ ’* eStab,1shma " t of acceptable strain 
points !„ ' tVI>eS - 8atWeen adJacent ‘train calculation 

would not Z\l t i and Pr ° babl> ° Ver tba antire • strain gradient 

ola ti ° Se " tb3t encompass an dastic-plastl c-elastic or a 

P astlc elastic-plastic variation. Therefore, 


if 



2 e yield I » 


remesh this element. 


Additionally, if the thermal strain changes In sign. Therefore, 
if 

e T 

”t < 0 
c 

j 

remesh this element. 
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Once the nonlinear solution has been entered, the element level checks 
become more complex and more Important. The total strain is now made up of 
the elastic strain, thermal strain, plastic strain, and creep strain: 


C 

i 



+ c 


c 

i 


Now stress and strain are no longer linearly connected; stress is a 
function of elastic strain only. Once again, between any two adjacent 
calculation points within one element, an elastic strain gradient greater than 
the allowable material elastic gradient is not desirable. Thus, 

if 


C . ~ C . ^ I | , 

i J 


remesh this element. The limit on the thermal strain would still be 
retained. 
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n 

If 


i 

\ ! 

e T 

4 - <° . 

e 

j 

remesh this element. 

j 

i 

The next check Is on the computed plastic and creep strain. No sign 
changes In either of these are allowed. add1t1on , . 9rad1ent „ 

1 

' 1 

if 

c p 

7 <0 

j 


n 

n 

•/ 

■ | 

or 



€ C 

-±< 0 
c c 

j 


i 

v ! 

1 

or 


, 1 

| e p - e p | > c 

i J P 


■ i 

1 

or 


| 

1 e c - c c | > C , 

i j c 


r ' — t 

** f 

1 



n 

remesh this element. 


i 

Next, proceed to the 
nature as the above, but 
elements. 

Interelement level check. These are of the same 
now Involve adjacent calculation points In adjacent 

' i 
\ 

f— V 

7 

' l 

"1 

i 
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2.8.3 Self-Adaptive Solution Strategies 


In the development of basic self-adaptive solution strategies, we have 
used the work of Edward T. Wilson of the University of California at Berkeley, 
and Joseph Padovan and Surapong Tovichakchaikul of the University of Akron. 

Wilson's efforts are directed toward an overall solution strategy, while 
Padovan's and Tovichakchaikul 's work is on load incrementing and time-stepping 
for geometrical and material nonlinear solutions. 

Wilson's philosophy on Internal program organization for SAP-80 computer 
programs Is applicable to the Component-Specific Modeling Program, with some 
extensions. He suggests that the basic Internal organization of a computer 
program for structural analysis depends strongly on the method used to form 
and solve linear equations, with the frontal and profile (or active column) 
methods most often used. Both have the exact same economy so that the choice 
must be based on other factors. 

In the frontal method, element stiffnesses and solutions of equations are 
formulated In a joint sequence manner. Therefore, all element stiffness 
subroutines, the equation solver, and the front of the stiffness matrix must 
be in core storage (or rolled in and out) during the reduction of the 
stiffness matrix. For the profile approach, the formation of all element 
stiffnesses for a particular type of element can be accomplished by a single 
call to one program segment. The formation of the total stiffness is a 
separate program segment in which the element stiffnesses are read In sequence 
from secondary storage and the total stiffness matrix is formed in active 
column blocks. In this case, the actual solution phase is another separate 
program link. Evaluation of substructure stiffnesses, calculations of mode 
shapes and frequencies, and evaluations of reactions and member forces are all 
separate links. This clear uncoupling of different phases of the program 
gives the profile approach a clear advantage in modularity and adaptive 
solution techniques. Also, the profile approach has no significant 
disadvantages when compared to the frontal method. 
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Padovan and his co-workers at the University of Akron have been developing 
Self-Adaptive Incremental Newton-Raphson Algorithms" for nonlinear problems. 
They use a three-level approach. In the first level, incremental Newton- 
Raphson operators are used to "tunnel" Into the problem solution space. The 
second level involves the constant monitoring of the different stages of 
solution through various quality/convergence/nonlinearity tests. The third 
level works with the results of the second level. The violation of any of the 
quality/convergences/nonlinearity tests triggers various scenarios for 
modifying the Incremental Newton-Raphson strategy. The self-adaptive 
modifications triggered by the third level fall into one of three categories: 
global stiffness reformation; preferential, local reformation; or load 
Increment adjustment. Recently, they have developed constrained, 
self-adaptive solution procedures for structures subject to high temperature 
elasti c/plastic/creep effects. In this, they used closed, piecewise, 
continuous least-upper-bounding constraint surfaces that control the size of 
successive dependent variable excursions arising out of the time-stepping 
process. 

A list of parameters to be controlled by the self-adaptive solution 
strategies has been generated. The parameters defined to date are listed 
below. 

Parameters to be Control led 

1. Element Type(s) 

2. Type(s) of Integration 

3. Order(s) of Numerical Quadrature 

4. Maximum Number of Iterations 

5. Tolerance(s) on Convergence 

6. Constitutive Equation(s) 

7. Yield Criterion (Criteria) 

8. Load Increments 

9. Time Increments 

10. Nonlinear Solution Algorithm(s) 



H 

I '■ 


Experience will, the m-house programs has given us a goon basis for 
developing the necessary tolerances on convergence. First, convergence Is 

ITT ,'°f ^ " 0t 9,0ba " yi ^ U ,s •»'«*«. at each element or 
should h"! tVr e9rat,0n POlnt ' SaC0 " d ’ f ° r n ” erUal conditioning, limits 
for tL d. WMCh 1neUSt,C Stra '" S Sre to be zero. Third, 

el;?! * CtS ’ b ° th t em P er ature and stress cutoffs should he 

h? th T T t1me 9ependent 1ne ' ast,c «"’» 15 considered to be 
following convergence criteria for Incremental analysis are the 


Time Independent 


If 


"j e P < PCUTOFF, c p = 0.0 


then 


Ac pr < TOL = CONVERGENCE 
or 


Aepi ~ Aep (x-!) 
Ac pi 


< TOL - CONVERGENCE 


Time Dependent 


If 


TEMP < TOLC, e c = 0.0 
and/or 

°e < V C c 5 °*° 
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or 


Ao gI < CTOL = CONVERGENCE 


*°el to .(l-1) 

Ao ei 


< CTOL ■ CONVERGENCE 


The different convergence criteria are dictated by the wide material strength 
levels encountered in nonlinear analysis. We have also discovered that it is 
advantageous to be able to change these criteria during the course of an 
incremental analysis. 


One approach taken in nonlinear computer codes is the right-hand-side 
technique, in which the plasticity is accounted for by adding an additional 
force vector to the right-hand side of the system of equations. 

[K] (d> = {F} + {f p } 

The basic logic is as follows: 

1. Solve for displacements from 

[K] {d} = {F} + {f p } 

2. Using the displacements and the constitutive equations, determine 
elastic and plastic strains for each element. 

3. Check convergence. 

4. Make an estimate of plastic strains that will satisfy the 
constitutive equations, equilibrium, and compatibility. 



n 


r" 




5. Based on the estimate of plastic strains from Step 4. form a new 
plastic load vector and go back to Step 1. 

This Iteration scheme continues until the convergence criteria are satisfied. 

The plastic Iteration accounts for a considerable portion of the total 

C0mputer c ° st * nonlinear finite element code. Subst nt a, 

'T b r ‘—.mg the convergence of plastic 

ar v , r 9 th6 " t,Mte ° f the S0,Ut, °" ** < -Tee options 


calculation ^f e ’* Slmp1est ' **. current 

of the sit „ I m' IT ^ C °" St1tUt,Ve as the estimate 

schemes. ‘ US “ a ' n,etf>0d on right-hand-side Iteration 


r~ 


The second scheme is a modification of 
is essentially a successive-over-relaxation 

solution is given by: 


the original iteration scheme, and 
( SOR) scheme. The estimate of the 


c 


P 



a 



* current estimate of solution 

“ previous calculate of plaatic strain from constitutive equations 

- current calculation of plastic strain from constitutive equations 

current acceleration factor 

= 1.5 is used 


r— 
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This estimation procedure continues until c 1 < c 1 " 1 then th- 
following Is used: p p * 


e 1 - + 0.5 cp - t 1 "' 




flMd SC Ir 15 b " ed °" an A,tkM ' S ».*1. for a 

it. ♦ « * terat 0n ' A,thou a fl we are not really doing a fixed-point 

teratlon, the finite element equations behave In much the same way. The 
equation used In estimating the solution Is- 


*1 

c 

p 


i-2 


(.‘- 1 - .‘- 2 ) 2 

- w p L 

e 1 ♦ 2c 1 " 1 ♦ « 1_2 
P P P 


1 

) 

n 

i 

t 

n 


Where the symbols are as before and = calculation of plastic strain 
om the constitutive equations two iterations ago. 

The A1 then’s extrapolation works best when performed every third 
Iteration, in between Altken's extrapolations 

_ i 

C ■ C 

P P 


n 
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2,8,4 COSMO SYSTEM FOR COMBUSTOR NUGGETS 


n 

i ' 

i 


r- - 
i 


r" 


Figure 33 shows a flow chart of the overall COSMO system Including the 
action positions of the adaptive controls developed In this program. For the 
combustor, the following adaptive controls have been Incorporated Into the 
system (the numbers are consistent with Figure 33 ). 


1. time increment 

2. load increment 

3. plasticity tolerances 

4. creep tolerances 

5. number of master region elements 

6. number of slices 

7. position of slices 

8. row refinement 

9. element refinement 


bnino a ! adapt1 ve contro ^ are a function of the structural coda 

un a b i ° r r! h,S SySte " the C0<)e an<l the COntro ' s •« «.«. developed 
ondar 30 Inalastlc Analysis Mathods for Hot Sactlon Structures.- Tha othar 

a aptlva controls ara keyed from a decision grid as Indicated In Figure 34 

The gradients In normalized stress, total strain, plastic strain, and creep 

strain will be used to rank requirements. 


r~- 


r~ 
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Figure 33. System Flow Chart Showing Adaptive Control Positions. 
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Figure 34. Combustor Nugget Decision Grid. 
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3 *° EXECUTIVE summary 


~~^=-~=r 

rr :,;r 

grated system for accurate and efficient prediction of temperature 
mlss^ 1 ’^ Stf T an<l Stra1 " h1stor1es thr oughout a complete flight ’ 

:::;:-va:;.r e ::;i;n;r d ,or ~ — ■«— 

as depIcted^rTFIgure^l "TcT ““ aCC0 " p,1shed 1n tw » Parallel work efforts 
existing .nalyslslthlds^d.^ o;";.;::^!? 0 ," MalUat, °" 

and development of geometric modeling and structure" 

rr, :=rrjrs:t :,rrr 

decomposition and synthesis capability. 
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APPENDIX A 

TASK II - DESIGN OF STRUCTURAL ANALYSIS 
SOFTWARE ARCHITECTURE 
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PREPROCESSOR ATTRIBUTES: 


A. SOFTWARE 

i 

• Modular in Structure 

n 

• Machine Independent 

i 

• Low Core Requirement 

* — i 

• Extensive Documentation 

j 

B. FEATURES/OPERATION 

n 

• User Friendly 

i 

• Extensive Diagnostics 


• Complete Data Summary 

! 

i 

• Interactive Graphics 

n 

. ) 

- Model Geometry 

- Material Physical Properties 

* — \ 

i 


n 

r— \ 
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n 
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POSTPROCESSOR ATTRIBUTES: 


A. ' SOFTWARE 

• Modular in Structure 

• Machine Independent 

• Low Core Requirement 

• Extensive Documentation 

B. FEATURES/OPERATION 

• User Friendly 

• Will Function in Batch or Time~Share Environment 

• Will Have Extensive Graphics Capabilities 

- X-Y Graph Plots 

- Contour and Deflected Shape Plots at Defined Planes 

• Will Have Built-in Data Manipulation Routines (User Controlled) 

- User Friendly Data Base 
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APPENDIX B 


COMPONENT TEMPERATURE AND PRESSURE DECOMPOSITION 
AND SYNTHESIS PLAN 


I • Blades and Vanes 


A. Input will be : 


1 . A series of mid~span station temperatures at a specified SS 
r *f er *oce case condition, and for each station a percent cord 
envelope dimension. 

2. A set of output percent cord envelope dimensions. 

3. A table of percent radial span (or radii) versus temperature 
factors, RF, where: 


RF 


Ts-T3 

Tms-T2 


Ts = temperature at specified span dim. 

Tms = temperature at midspan. 

The point density in this table will be such that spanwise linear interpo 
lation will suffice. 


B. Calculations 

1. For each temperature a nominal cooling effectiveness will be 
calculated from: 

_ T4l-Tm 
n c “ T41-T3 


2. At each new condition, all station r| values will be modified 
by a factor, CF, as follows: C 


CF = 


( l "H c ) 

(l-ri c )REF 


= (T3'*T4l) 0f /(T3*T4l) 


a 

REF 


where a is input 


3. Output Nc values will be linearly interpolated based on percent 
cord envelope at adjacent input stations. 

4. For each output radial distance (specified in the input) all 
station temperatures will be calculated as follows: 
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Tij • T3 + (l-n CR )*CF(T4;-T3)*RF 


Where i = station index 

j * radial station index 


II. Combustor 


A. Inputs will be ; 

1. Metal bulk made temperatures at axial stations on inner and 
outer liners, and for each station, x and y dimensions of 
node centers at a specified reference case condition. Both 
hot streak and average metal temperature values will be input. 

2. It is assumed that output node dimensions will match the input 
nodes. If cross-meshing is to be required, it will be done 

before the input is defined. 

• 

3* A set of linear equation constants will be input for a sparse 
set of* axial locations, identified by AT station: numbers. 

Output AT locations will be same as input locations. 

4. A set of AP scaling constants will be input at a sparse set of 
locations, identified by AP station numbers. Output AP sta- 
tions will be same as input. 

B. Calculations 


1. A nominal cooling effectiveness will be calculated for each input 
combustor average and hot streak temperature. 

2. At each flight condition, each metal temperature will be recalcu- 
lated from: 

T “ij = T3 j + d-V ) (Ul.-H.) 

Where i * station index 

j ■ flight phase index 


3. At each flight condition, each AT will be recalculated from: 
AT ij = C(l-n c )*T41 + (n c -1)*T3 ] [b +m *P3.] 

^ J -Llj 

4. At each flight condition, each AP will be recalculated from: 

AP = P3 *K * 2 * T3 

ij j i P3 j j 
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APPENDIX C 

THERMODYNAMIC ENGINE MODEL SPECIFICATION 
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I 


Input 


II. 


III. 


A. Setup information (furnished with model) 

1. Engine performance data; 148 cases per Table I-A-l. 

a. Case parameters per Table I-A-l-a. 

2. Engine rating data - Table I-A-2.1-3. 

3. Power level index matrix, Table I-A-3. 

B. User Information 

1. Mission definition data 

a. One line for each mission phase point. 

b. Each line contains the following control data: 

1. Phase # 

2 . Mach number 

3. Altitude - feet 

4. Offset from standard day - *F 

5. Power level parameter code # 

6. Power level parameter value 

7. Customer bleed - #/sec. 

8. Deterioration level - *F 

9. Time increment between this phase point and the 
next, min. 

c. One line, following the mission phase point data line, 
for each parameter to be offset from its steady state 
value . 

d. Each offset line shall contain: 

1. Phase # 

2. Parameter # 

3. Offset factor 

4. Offset adder 

Output 

A. A performance case for each mission phase point 

1. Parameters per Table I-A-l-a. 

2. Format similar to Table II-A-2. 

Technical Basis 


A. Each new case will be generated from available cases (I-A-l) by 
a disciplined interpolation process similar to that currently 
used in the Life Analysis by Stress and temperature Simulation 
(LASTS) program. 

1. All parameters will be transformed to a functional form that 
has optimal linearity relative to all other parameters, 
a. A study will be performed on CF6-50C2 engine 
performance data to evaluate and improve the 
interpolation functions. 
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2. Each transformed parameter will be interpolated by a 
linear interpolation process and then transformed back 
to its normal fora. 

3. The interpolation "targets" shall be specified in the 
input for each mission phase point (I-B-l-b). 

4. The interpolation proeess for each phase point shall begin 
with a base case n«ar the desired mission phase point. 

5. Each interpolation step will convert the base ease (or 
previously modified case) to the desired level of the 
target control parameter (i.e., MN) . 

a. Linear partials are assumed; interactions are ignored. 

b. Partials will be derived from two or more "partial 
cases" near the base case conditions. ' 

c. The interpolation steps will be performed in the 
sequence in Table I-A-l-a. 

d. The specific power level parameter to.be used as an 
interpolation target is input as "power level 
followed by the target value. 

e. For flight idle and ground idle, a special power 
level # (one for FI or one for GI) will be entered, 
followed by a 2ero parameter value; the standard 

FI and GI power levels will be used. 

f. For thrust reverse, a special power level # will be 
used and a value of fan speed vill follow; thus 
thrust reverse power level will always be based on a 
fan speed target, 

The LASTS interpolation process will be modified to eliminate 
the current manual procedure for the generation of interpolation 
instructions. 

1. A set of base case numbers will be provided as a function of 
altitude and Mach number. * 

2. A family of pairs or triplets of partials case numbers will 
be provided as a function of altitude and Mach number for 
each control parameter. 

3. The user will be required to input only the data in I-B-l-b 
for each mission phase point. 

The Thermodynamic Engine Model (TOE) shall have the capability 
to predict the minimum time for speed changes due to throttle 
actions . 

1. The user shall have the option to input zero throttle-action 
transient times, and the model will calculate appropriate 
transient times, subtracting them from following phase 
times. 

2. The transient time calculation shall be sensitive to the 
effects of altitude. 

The user shall have the option of selecting an appropriate 
CF6-50C2 power management point, and avoid the need for 
specifying absolute values of the power level parameter. 

1. The power management parameter code shall call for take-cff, 
max climb, or max cruise rating, and the adjacent value 
shall specify the I derate desired based on Z thrust. 


E. Offsets of specific parameters, relative to the steady state 
performance cases shall be permitted to simulate take-off 
transient conditions. 

1. Each parameter change shall be specified by a line 
following the mission phase point data. 

2. The parameter offset data shall be: 


a. 

Case # 


b. 

Parameter # 

c. 

Offset 

factor 

d. 

Offset 

adder 


3. Offset calculations will be performed after mission phase 
point interpolations are completed. 

Software Characteristics/Interfaces 

A. Later (to be integrated with overall software of the COSMOS 
Program) . 
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TABLE I-A-l 

CF6-50C2 PERFORMANCE CASES 


M 


A PCNLR PCN5R2 
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ATo WB DT49 °F 
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TABLE I-A-la 

ENGINE PERFORMANCE CASE PARAMETERS 


P2 

Fas Inlet Total Pressure 

PSIA 

P3 

Compressor Discharge Total Pressure 

PSIA 

P4 

Turbine Inlet Total Pressure 

PSIA 

P49 

Turbine Outlet Total Pressure 

PSIA 

T2 

Fan Inlet Total Temperature 

*F 

T3 

Compressor Discharge Total Temperature 

*F 

T41 

Turbine Inlet Total Temperature 

*F 

T49 

Turbine Outlet Total Temperature 

* F 

V25 

Fan Air Flow 

#/sec 

FNIN1 

Installed Thrust 

# 

DTAMB 

Offset from Standard Day Temperature 

•f 

W41 

Turbine Air Flow 

#/ sec 

XNH 

Core Speed 

8PM 

XNL 

Fan Speed 

RPM 

MN 

Mach Number 


ALT 

Altitude 

Feet 

WB27/VB3 

Customer Bleed 

#/sec 

DT49 

Engine Deterioration Index 

# F 
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TABLE I-A-2.1 

APPROX. CF6-50C2 RATING DATA 
TAKEOFF 


°C Feet 

T2 ALT PCNLR XNLR ALT 


-60 0 107.75 3736 1000 

20 
23 

25.5 
27 

20.5 

30.5 


50 

105.7 

3620 

60 

97.8 

3357 


T2 

ALT 

PCNLR 

XNLR 

ALT 

-60 

3000 

111.85 

3039 

4000 

20 





23 





25.5 





27 


111.1 

3814 


20.5 


110.4 

3789 


30.5 


109.75 

3767 


50 


105.7 

3620 


60 


97.8 

3357 



PCNLR 

XNLR 

ALT 

PCNLR 

XNLR 

110.4 

3789 

2000 

111.1 

3814 




110.4 

3709 

109.75 

3767 

109.75 

3767 

105.7 

3628 

105.7 

3628 

97.8 

3357 

97.8 

3357 


PCNLR 

XNLR 

ALT 

PCNLR 

XNLR 

112.75 

3070 

* 5000 

113.82 

3907 




112.75 

3070 

111.05 

3839 


111.85 

3039 

111.1 

3014 


111.1 

3814 

110.4 

3789 


110.4 

3789 

109.75 

3767 


109.75 

3767 

105.7 

3620 


105.7 

3628 

97.8 

3357 


97.8 

3357 


J 


: ] 




] 


j 


j 


j 



°c 

FEET 



T2 

ALT 

PCNLR 

XNLR 

-60.0 

-7.0 

36009 

117.1 

4019 

-5.7 


116.2 

3989 

-1.2 


113.6 

3099 

11.0 


106.7 

3662 

21.2 


103.4 

3549 

34.0 


90.5 

3301 

60.0 


90.4 

3103 


11 

ALT 

PCNLR 

XNLR 

-60.0 

-7.0 

-5.7 

20000 

106.7 

3662 

-1.2 




11.0 




21.2 


103.4 

3549 

34.0 


98.5 

3301 

60.0 


90.4 

3103 


1 


TABLE I-A-2.2 

APPROX. CF6-50C2 RATING DATA 
MAX CLIMB 


FEET 

ALT 

PCNLR 

XNLR 

42000 

116.2 

3989 


113.6 

3899 

106.7 

3662 

103.4 

3549 

98.5 

3381 

90.4 

3103 


ALT PCNLR XNLR 

IOOOO 103.4 3549 


98.5 3301 

90.4 3103 


FEET 

ALT 

PCNLR 

XNLR 

30000 

113.6 

3899 


106.7 

103.4 

98.5 

90.4 

3662 

3549 

3381 

3103 

ALT 

PCNLR 

XNLR 


98.5 

3381 


90.4 

3103 


TABLE I-A-2.3 

APPROX. CF6-50C2 RATING DATA 


to 

00 


MAX CRUISE 


°c 

FEET 



FEET 

T2 

ALT 

PCNLR 

XNLR 

ALT 

-4.5 

36089 

112.42 


30000 

- 1.0 


109.70 



0 


109.25 



9.2 


104.45 



♦20 


100.64 



♦29.2 


97.47 



♦43.7 


92.70 



♦60.0 


B7.35 




T2 

ALT 

PCNLR 

XNLR 

ALT 

-4.5 

25000 

104.45 


20000 

-1.0 




0 





9.2 





20.0 


100.64 



29.2 


97.47 



43.7 


92.70 



60.0 


07.35 



T2 

ALT 

PCNLR 

XNLR 


-4.5 

0 

92.70 



-1.0 





0 





9.2 





20.0 





29.2 





43.7 





60.0 


07.35 




PCNLR 

XNLR 

FEET 

ALT 

PCNLR 

XNLR 

109.70 


42000 

109.25 


109.25 

104.45 

100.64 

97.47 

92.70 

87.35 



104.45 

100.64 

97.47 

92.70 

87.35 


PCNLR 

XNLR 

ALT 

PCNLR 

XNLR 

100.64 


10000 

97.47 


97.47 

92.70 

87.35 



92.70 

87.35 



J __] 1 J .... 1 _ J J 
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TABLE I-A-3 


AUTOMATED INTERPOLATION SYSTEM 
POWER LEVEL INDEX 
CF6-50C2 


P.L. INDEX 

m — 

O 

PCNLR 

ALT 

10000 

1 

109.0 


119.0 

2 

100.0 


109.0 

3 

90.0 


98.0 

4 

75.0 


81.0 

5 

55.0 


58.0 

6 

38.0 


44.0 

7 

24.0 


34.0 

P.L. INDEX 

ALT 

10000’ 


ALT 

40000’ 

1 

103.5 


120.0 

2 

91.5 


106.0 

3 

75.0 


85.0 

4 

58.0 


67.0 

5 

44.0 


48.0 

6 

34.0 


37.0 
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TABLE II-A-2 

PERFORMANCE OATA OUTPUT FORMAT 


" ENGINE PERFORMANCE OAT* BY MISSION PHASE • 


PHASE CASE ** ENGINE PARAMETER • • 

0 0 P4 * PS was FNIN1 


0TAM8 


1 197 

2 190 

3 1*1 

4 1S2 

9 193 

9 191 

7 191 

• 192 

* 193 

10 154 

11 19S 

12 159 

13 197 

14 187 

19 188 

16 189 

17 180 

18 19 1 

19 192 

20 193 

21 194 

22 194 

23 165 

24 189 

29 187 

29 187 

27 194 

28 195 

29 198 

30 186 


14.478 14.478 1.108 Q 

14.476 14.476 1 . 108 Ol 

14.476 14.478 1.108 O. 

19.178 14.476 1.108 0. 

-18,476 18.476 1. 1Q8 q, 

16.912 14.895 38.502 1943.746 

16.912 14.99S 38.502 1943.746 

80.731 20.930 243.413 42726.324 

82.839 21.907 245.733 39699.600 

66.172 18.773 206.290 30540.937 

S3* 4 * 7 16.309 194.796 21481.488 

34.321 8.196 109.120 9766.940 

,24.340 9.9 19 81.158 6421.298 

24.3*0 5.915 . 81.188 9421.258 

4.888 2.797 19.617 *892.989 

6.888 26.739 -'1914.288 

,13.099 10.748 39.909 *2310.870 

16.297 13.944 38.100 894.115 

27.946 14.707 95.930 7275.440 

29.131 15.402 99.315 8873.031 

22.015 14 , 86 1 68.412 4493.193 

22.015 1'4.86t 68.4 >2 *493.193 

43.788 16.292 145.437 14871.125 

43.084 16.184 144.799 19603.446 

16.9*2 14.622 38.588 19*6,327 

16.942 14.622 38.958 19*6.327 

16.503 16.903 0.974 0. 

14.903 14.503 0.974 0. 

14.903 14,503 0.974 p, 

14.503 14.503 0.974 0. 


1.7QQ 

1.700 

1.700 

1.700 

1.7QQ 

1.700 

1.700 

1.701 

2.501 

2.502 
9.001 

12.200 

12.200 

12.200 

13.200 

7.800 

4.200 

2.500 

2.500 

1.700 

1.700 
! . 700 
1.704 
1.700 

1.7QQ 

1 . 700 
1.700 
1.700 
1 ■ 700 
1.700 
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